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ABSTRACT 

We explore the variability and cross-frequency correlation of the flux density 
and polarization of the blazar OJ287, using imaging at 43 GHz with the Very 
Long Baseline Array, as well as optical and near-infrared polarimetry. The polar- 
ization and flux density in both the optical waveband and the 43 GHz compact 
core increased by a small amount in late 2005, and increased significantly along 
with the near-IR polarization and flux density over the course of 10 days in early 
2006. Furthermore, the values of the electric vector position angle (EVPA) at 
the three wavebands are similar. At 43 GHz, the EVPA of the blazar core is 
perpendicular to the flow of the jet, while the EVPAs of emerging superluminal 
knots are aligned parallel to the jet axis. The core polarization is that expected 
if shear aligns the magnetic field at the boundary between flows of disparate 
velocities within the jet. Using variations in flux density, percentage polariza- 
tion, and EVPA, we model the inner jet as a spine-sheath system. The model 
jet contains a turbulent spine of half- width 1.2° and maximum Lorentz factor 
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of 16.5, a turbulent sheath with Lorentz factor of 5, and a boundary region of 
sheared field between the spine and sheath. Transverse shocks propagating along 
the fast, turbulent spine can explain the superluminal knots. The observed flux 
density and polarization variations are then compatible with changes in the di- 
rection of the inner jet caused by a temporary change in the position of the core 
if the spine contains wiggles owing to an instability. In addition, we can explain 
a stable offset of optical and near-IR percentage polarization by a steepening of 
spectral index with frequency, as supported by the data. 

Subject headings: polarization — BL Lacertae objects: general — BL Lacertae 
objects: individual (OJ287) 



1. Introduction 



Blazars are powerful and variable sources of polarized radiation across a wide range 
of wavelengths. However, establishing the location for this energetic and often dynamic 
emission within blazars is problematic. In particular, efforts to determine the site of op- 
tical emission within the parsec-scale jet have been inhibited by the low spatial resolution 
available at this waveband. Recently, new techniques have been developed to counter these 
difficulties. Unresolved optical emission can be connected to a location in resolved Very 
Long Baseline Array (VLBA) images of the jet by comparing electric vector position an- 
gles (EVPA) and degrees of polarization, as well as variable total flux density across the 
two wavebands. Studies utiliz ing this method have concluded that variable optical emission 



originates in the 43 GHz c ore (ILister fc Smithll2000l : iGabuzda et al.ll2006l ; iD'Arcangelo et al. 



20071 : IJorstad et al.l 120071 ). Here, we extend this technique through time- variable correla- 
tions between optical and radio wavelengths in the blazar OJ287 (0851+202, z = 0.306; 
Sitko &: Junkkarinenlll985l ). 



O J287 is a BL Lacer tae object contai ning a one-sided jet with complex structure on par- 
sec ( [Roberts et al.lll987l ) and kiloparsec (IPerlman &: Stockd Il994l ) scales. Previous studies 
have imaged the jet of OJ287 at milliarcsecond (mas) resolution at multiple radio wave- 
lengths; the resultant maps contain a compact core with substantial linear polarization 
(P< 10%), and a tendency toward a 90° difference between the EVPA of the core and 



the EVPA of th e more extended jet ([Roberts et al.l 119871 ; IGabuzda fc Cawthornd Il996 
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position angle as cyclic, with period 11.6 years. Such periodic variation, plus projection 
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effects, agree with the tendency of OJ287 to establish moderately long-term trends, such as 
the observed position a ngle of ~ —110° in the late 1 980's and ~ ~90° m t he mid 1990's 



(Kol 



2001 



2005). In addition. 



gaard et al.l Il992l ; iGabuzda fc Cawthornd Il996l ; iTateyama et al.l Il999l ; IJorstad et al. 



Jorstad et al.l (120051 ) noted the presence of several stationary knots 
in the jet of OJ287 with position angles ranging from —90° to —111°, a possible lingering 
effect of the swinging inner jet. 

The inclusion of data at optical waveleng ths has established a preliminary connectio n 
between the o ptical and radio EVP A in OJ287 dfludnick et al. 1978 ; Gabuzda fc Sitko 1994 ). 
Furthermore, iKikuchi et al.l (119881 ) and ISillanpaal (Il99ll ) have measured rotations of EVPA 
at both wavebands, over scales varying from a month to a decade. 

We expand previous studies of OJ287 by combining resolved 43 GHz VLBA images 
of the parsec-scale jet with simultaneous polarimetric measurements at optical and near- 
IR frequencies over two intensive 10-day campaigns. Our results show a clear connection 
between the polarized emission at optical and near-IR wavelengths and that associated with 
the 43 GHz parsec-scale core. In addition, we find rapid, synchronous changes in radio, 
optical, and near-IR degree of polarization and total flux in early 2006, which we interpret 
as the consequence of a changing inner jet relative to the line of sight. 



2. Observations 

In our observational program, we utilized instruments at optical, near-infrared, and radio 
wavelengths to collect highly-sampled measurements of percentage polarization, EVPA, and 
flux density. We obtained this information for 21 blazars during coordinated campaigns in 
2005 October-November and 2006 March- April. 

To maximize the variability sampling, we observed at optical wavelengths using two 
facilities. The majority of the data are from the 1.55 m Kuiper telescope at Steward Ob- 



servatory, equipped with the SPOL spectropolarimeter (ISchmidt et al.lll992l ). We made ten 
spectropolarimetric observations of OJ287 between 2005 October 25 and November 3 and a 
further eight observations between 2006 March 27 and April 4. For these 18 measurements, 
we utilized either a 3.0" or 4.1" slit width and a 10" spectral extraction aperture. A diffrac- 
tion grating with 600 lines per mm was used, producing a diffraction order of 400 to 800 
nm with a dispersion scale of 0.4 nm pixel -1 . We calculate the spectral resolution to be 
approximately 1.9 nm for a 3.0" slit width and 2.5 nm for a 4.1" slit width. We binned the 
spectropolarimetry to derive polarization measurements within wavebands similar to those 
of the B and V Johnson filters and R and / Kron-Cousins filters. This binning represents 
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the continuum, as no line features were present in the spectra. Photometric conditions ex- 
isted on nine nights in 2005 and on two nights in 2006. On these occasions, we measured 
the optical flux using either a 5.1 x 12", 7.6 x 12", or 7.6 x 14" extraction aperture during 
the 2005 campaign. We employed a 7.6 x 12" extrac tion aperture in 2 006. We calibrated 
by comparison with field stars using magnitudes from ISmith et al.l (119851 ). A more detailed 



descri ption of the procedural aspects of the spectropolarimetry can be found in ISmith et al. 
f)2003h . 



We also performed optical observations at the Crimean Astrophysical Observatory on 
five nights between 2006 March 26 and 2006 April 5. We took measurements on the AZT-8 70 
cm telescope equipped with an ST7-based photometer-polarimeter, in which Savart plates are 
aligned at 45° to one another. The orientation of the plates thus produces a measurement 
of either the Stokes Q or U parameter. Up to five images in each position of the plates 
allowed co-addition to improve the signal-to-noise ratio. We obtained both polarimetric and 
photometric data in the R waveband on each night using aperture photometry, and calibrated 
using standard field stars. We also include in this study monitoring at this observatory 
outside the 2006 campaign period. OJ287 was not observed from the Crimean Astrophysical 
Obse rvatory during the 2005 cam paign. All optical measurements are corrected for statistical 
bias fjWardle fc Kronberglll974l ). 



We observed in the near-infrared H-ban d (1.65 iim) at Lowe ll Observatory's Perkins 1.8 
m telescope. Using the Mimir instrument (jClemens et al.l 120071 ). we obtained polarimetric 
and photometric measurements twice during the 2005 campaign and on 3 nights during 2006. 
Polarimetry involved imaging OJ287 through a rotating half-wave plate and a wire-grid an- 
alyzer, with 32 images obtained over a full 360° in position angle of the wave plate. Images 
were calibrated using a bad pixel mask, dark current and bias count subtraction, flat-field 
correction, a correction for the non-linear response of the detector, and EVPA calibration 
based on observation of polarimetric standards. In addition, we performed differential pho- 
tometry relative to nearby IR standard stars. 

We observed OJ287 with the VLBA at 43 GHz on three dates spaced approximately 
equally across the 10 days of each campaign: 2005 October 24, October 28, and November 
2, and 2006 March 27, March 31, and April 5. Correlation of the data collected occurred 
at the Array Operations Center of the National Radio Astronomy Observatory (NRAO) in 
Socorro, New Mexico. We calibrated the uv data using the Astronomical Image Processing 
System (AIPS) software from N RAO, and then constr ucted and ana l yzed the images with 
the DIFMAP software package (jShepherdl 119971 ). See IJorstad et al.l (120051 ) for a complete 
discussion of our data reduction methods. The EVPAs in the images were calibrated by sum- 
ming the total and polarized emission across images of a number of objects, and comparing 
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these integrated polarization characteristics with those supplied by NRAO from measure- 
ments made with the Very Large Array (VLA) at dates as close as possible to our epochs 
of observation. The objects used, with the average separation in time in parentheses, were 
0923+392 (4C +39.25) (6 d), BL Lac (3 d), and OJ287 (6 d) for the 2005 campaign, and 3C 
273 (7 d), 3C 279 (7 d), OJ287 (7 d), and B2 1156+295 (4C +29.45) (7 d) during 2006. The 
variance in the EVPA calibration, as determined using the integrated VLA data, was used to 
estimate the uncertainty in the calibration. To check the validity of our EVPA calibration, 
we examine the change in instrumental D-terms between epochs, as well as the variance in 
EVPA of jet components with known values. We conclude from this examination that the 
EVPA calibration is correct to within ±5° in 2005 and to within ±6° in 2006. We apply a 30° 



compe nsation for Faraday rotation to the core region of O J287, as derived by iJorstad et al. 



(120071 ); unless otherwise noted, all EVPAs quoted have been corrected for this Faraday rota- 
tion. To determine the motion of components within the jet, we use supplementary 43 GHz 
images of OJ287 observed with the VLBA at four epochs - 2006 May 16, June 17, August 
10, and October 5. These supplementary images are part of a monthly monitoring program 
of a sample of blazars, and have wf-coverage comparable to that obtained during the 2005 
and 2006 campaigns. 



3. Analysis 

The jet of O J287 consists of both superluminally and subluminally moving components 
that extend from the core along a position angle (PA) that varies among the components. We 
derive the total flux density, polarized flux density, position, and size of these components by 
fitting each /, Q, and U map with circular Gaussian components using the task MODELFIT 
in Difmap. Components are fit in order of peak flux densit y, from largest to sm allest. A 



detailed description of the fitting technique can be found in IJorstad et al.l (120051 ). Table 1 
lists the fitted parameters for each component - position relative to the core, PA, total flux 
density, degree of polarization, and EVPA - as well as the reduced chi-square statistic for the 
cumulative model after adding each component, starting with the core. In Figure [U we show 
a VLBA image for each of the two campaigns, as well as a VLBA image from each of the four 
supplementary epochs, with positions of three components marked. In the right panel, we 
utilize a small circular Gaussian restoring beam of 0.1 mas FWHM diameter, corresponding 
to the resolution of the longest baselines, to demonstrate the presence of components close 
to the core. 



In October 2005, the submilliarcsecond structure of OJ287 included the emergence of a 
new feature. This component is labeled in Table 1 and all figures as the 'West Component" . 
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In Figure [2] (top), we display the three high-resolution images of this component from the 
2005 campaign. These images reveal the presence of a strong core with flux density 2.2 ±0.2 
Jy, EVPA x — 19° ± 4° (after correction for Faraday rotation), and degree of polarization 
P = 12.1 ±0.5%, and a component with flux density 0.8 ±0.2 Jy to the west (PA -74 ±14°) 
with x — ~ 79 ± 10° and P = 6.9 ± 0.8%. We do not detect this component during the next 
campaign (2006 March). 

During the 2006 campaign, a component (labeled "2" in Figures [TJ [2J and [3]) separated 
from the core to the west at PA = — 97 ± 11°, with long-term motion along PA = —114 ±4°. 
We first imaged this knot during the 2006 campaign, and subsequently in 2006 May, June, 
and August. The emergence of the component in the three epochs of observation of the 2006 
campaign is shown in high resolution in Figure [2] (bottom). During the campaign, the new 
knot was positioned 0.1 mas to the west of the core, with a flux density of 1.8 ± 0.5 Jy and 
X = 77 ±4°, while the core had x = —1 ±4°, after correction for Faraday rotation. The knot 
retained a steady value of x — 66 ± 11° between March 31 and August 10. It s flux density 



dropp ed steeply and steadily, as found for three previous knots in O J287 by IJorstad et al. 



(120051 ). to a value of 0.14 Jy on 2006 August 10. 



A plot of separation from the core versus time (Figure [3l left) reveals that the proper 
motion of component 2 is approximately constant throughout its lifetime at 0.8 ±0.1 mas 
yr _1 . This corresponds to an apparent speed of 12±lc, and a "birth" (coincidence with the 
core) on 2006 February 4 ± 27 days0 

In late 2006, we observed a component moving along a more westward long-term path, 
-88 ± 14° (labeled "3" in Figures El and [3]). From the 43 GHz images of 2006 August 10 
and October 5, we derive an angular velocity of 0.7±0.2 mas yr _1 , or an apparent speed of 
10±3c. We then calculate a birth date for this feature of 2006 June 17 ± 33 days. Within 
the two dates of observation, the flux density of this component decreases by a substantial 
amount, as can be seen in Table 1. 

In addition, we map the path of an older, slower component along a position angle of 
— 147 ± 5° (labeled "1" in Figures [IJ [2J and [3]). Component 1 appeared in the images of 
both the 2005 and 2006 campaigns, as well as in 2006 May, June, August, and October. We 
observe it to have x — — 18±15° throughout its lifetime. Intermittent excursions in x , as well 
as substantial and significant excursions in P and flux density for this component suggest 
interaction with the surrounding medium. We calculate an angular velocity of 0.37±0.04 
mas yr _1 , or an apparent speed of 5.3±0.6c, which leads to a birth date of 2004 December 



1 We use the cosmological parameters H = 72 km/s/Mpc, Q, m = 0.27, and Qa = 0.73 to calculate all 
speeds and length scales. 
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13 ± 77 days. Figure [3] (right) shows the paths of the components born in 2004 December, 
2006 February, and 2006 June. 

Figure H] shows both ^(optical) and x(43 GHz core) during the 2005 campaign. We find a 
short rotation (from 20° to 8°), then plateau, of x(optical). The average x(-R-band) is 12±5°, 
the average x{H-band) is 8 ± 6°, and the average x(43 GHz core) is 19 ± 4°, comparable 
values. Within the 2006 campaign (JD 2453820 - JD 2453832), shown in Figured we observe 
a steady rotation of —2° day -1 in both x(optical) and x(H-band), with average values of 
x(-R-band)= ± 7° and x(iJ-band)= —1 ± 9°, while the 43 GHz core value is constant at 
x(43 GHz core)= —1 ±4°. Rotation in the latter is not apparent, perhaps due to the rather 
high uncertainties, but x(43 GHz core) agrees with the average x(°pti ca l) an d x(H-band) 
within the errors throughout the period of observation. 

The degrees of polarization in the optical and radio core are nearly steady throughout 
the 2005 campaign, as shown in Figure HJ with P(R) increasing by only ~ 0.1% day -1 . We 
measure the mean radio core polarization to be 12.1 ± 0.5%, the mean if-band value to be 
26 ± 2%, and the mean i?-band value to be 30.7 ± 0.8%, with a consistent ratio P(R)/P(A3 
GHz core) = 2.5. We also note that the optical polarization appears to oscillate with a period 
of roughly 2 days that is observed at all optical wavelengths (see Figure [6]). We observed 
stability in 5 polarized standard stars during these observations, as well as polarization 
< 0.1% for 2 unpolarized standards, leading us to conclude that this oscillation is not a 
systematic, instrumental effect. In addition, we not e that this behavio r is not observed in 



any other blazar during the campaign. Furthermore, ISmith et al.l (119871 ) find fluctuations in 
the degree of polarization of OJ287 of a similar period. We conclude, due to the repeating 
nature of this behavior, that the effect is intrinsic to the blazar, and not an artifact of 
propagation through intervening media. However, we do not have 43 GHz sampling on a 
similar timescale, so that the data are insufficient to draw any conclusions as to the structural 
basis for this variation. 

The flux density underwent a modest change during the 2005 campaign. There is a 
36% increase in i?-band flux density from 12.3 to 16.7 mJy and a less dramatic increase 
in 43 GHz core flux density from 2.0 to 2.3 Jy. We fit the o ptical spectra r esulti ng from 



the spectropolarimetry, corrected for Galactic extinction using ISchlegel et al.l ([1998), with a 
power law and find an optical spectral index, steady over the campaign, of a op t = — 1.28 ± 
0.05, where the total flux density F v oc u a . 

During the 2006 campaign, the optical, infrared, and radio core polarizations increased 
continuously (see Figure [5]), with no significant time delay between frequencies. The optical 
.R-band polarization began at 8.5% on 2006 March 27 and grew to 29.5% by 2006 April 
5. A least-squares fit to P versus time gives a growth rate of 2.30 ± 0.01% day -1 . The 
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infrared polarization increased from 10% to 22%, and P(43 GHz core) from 9.7% to 14.4%. 
We note that polarization at all wavelengths is approximately the same near 2006 March 
27, ~ 10%. After that date, the polarization increases at all wavebands, but the increase 
is greater at higher frequencies (see Figure [7]). The dominance of wavelength-dependent 
polarization effects grows with the percentage polarization. 

In addition, the optical, near-IR, and 43 GHz flux densities increased during the 2006 
campaign, as shown in Figure [5J In i?-band, we observed a steady doubling in flux density 
from 5.2 mJy on March 26 to 10.6 mJy on April 3 before it declined on April 4 and 5. The 
infrared if-band emission increased from 19.1 mJy on March 26 to 27.1 mJy on April 5. In 
the 43 GHz core, the flux density rose in a similar fashion, from 1.55 Jy on March 27 to 1.82 
Jy on April 5. From the optical measurements, which had a higher frequency of sampling 
than either those at 43 GHz or if -band, we conclude that the local maximum occurred on 
April 3 (JD 2453829). We also observed the optical spectral index to have a stable value of 
otopt = — 1.3 ± 0.1, which agrees within the measurement uncertainty with the index derived 
from the 2005 data. Combining the 5 if -band flux measurements over both campaigns, we 
calculate an average spectral index between if-band and J-band of olh-i — —0.9 ± 0.1. 



4. Discussion 

4.1. Multiwaveband Polarization Characteristics 

Since the dramatic increase in the degree of polarization of OJ287 during the 2006 
campaign occurs both at optical/near- infrared wavelengths and in the emission from the 43 
GHz core, it indicates that the regions of emission are at least partially cospatial. Further 
evidence for this is provided by the close correspondence of x(optical/IR) and x( r adio core) 
during both the 2005 and 2006 campaigns. We conclude that the variable optical/near- 
infrared emission originates from the core region of the parsec-scale jet. 

The polarization vector in the core of OJ287 is oriented roughly perpendicular to the 
flow of the jet. Under the standard relativistic synchrotron model for production of optically 
thin polarized emission - wherein \ is transverse to the magnetic field lines - the orientation 
of x i n OJ287 reveals that the mean magnetic field in the core is aligned parallel to the 
flow of the jet. To more precisely map the orientation of the magnetic field in the jet, we 
determine the closest observed epoch to the ejection of each component, and note x(core) at 
each of these epochs. We find a consistent difference of 99 ± 13° between x(core) at the time 
of ejection and the subsequent path of the emerging components. Therefore, the observed 
magnetic field in the core is approximately aligned with the axis of the jet at any given 
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epoch, parallel to the direction of flow. A region of velocity shear in the core is the simplest 
explanation for this observed field geometry. The shear could be produced when fast plasma 
near the jet axis flows past slower material closer to the boundary, resulting in a gradient 
of velocities that stretches and aligns the magnetic field, which may have originally been 
randomly oriented in turbulent cells, along the direction of flow. 

We also investigate the possibility that multiple radio components contribute to the 
change in polarization, and conclude that this is not a feasible solution. We observe a 
distinct stability of the orthogonality between x(core) and the jet axis, even for significant 
changes in the projected PA of the jet. Combined with the constant agreement between 
x(core) and x(R)> we conclude that the changes in \ and P are consistent with intrinsic 
changes in the viewing angle of the jet, which causes changes in the observed angle of aligned 
magnetic field. In addition, throughout the period of increasing polarization, we observe no 
variation in optical spectral index with time or frequency, nor do we observe any evidence for 
wavelength-dependent x- These are tw o tendencies that ar e commonly observ ed in the case 



of mu ltiple-component flares in OJ287 (IHolmes et al.lll984j ) and other blazars (IBrindle et al. 



19861 ). 



4.2. Model for Polarization Behavior 



During the 2006 observations of OJ287, we see the optical flux density increase, then 
decrease, while the percentage polarization continues to grow (Figure [5]). To explain this 
behavior, we create a structural model for the inner jet as a combination of three elements. 
The first element is a large, relatively slow turbulent sheath, which surrounds a thinner, 
faster, turbulent spine along the axis of the inner jet. The transition region between the 
two zones produces the final element, with the magnetic field there aligned longitudinally by 
shear. We use T = 16.5 as the average Lorentz factor of the spine, deri yed from observation s 
of the fast jet components of OJ287 in the detailed VLBA study of IJorstad et al.l (120051 ). 
The Lorentz factor of the bulk jet flow is calculated from the observed pattern speeds of 
components - which cannot exceed T - and the concurrent jet viewing angle, and is therefore 
greater than or equal to the fastest observed component speed. In addition, we set the 
half-opening angle o f the spine to be within the range 0.8 ± 0.4°, also derived for OJ287 in 
Jorstad et al.l (120051 ). The average Lorentz factor of the sheath is stepped through a series 
of values to find the best fit, and the boundary region is assigned a Lorentz factor that is the 
arithmetic mean between r spine and r s h ea th- As the sh eath is a large struc ture, we hold its 
viewing angle steady, using the value of 3.2° derived in I Jorstad et al.l (120051 ) . We then allow 
the viewing angle of the small spine and boundary region to vary with time. As the viewing 
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angle of a spine-boundary structure changes when the jet direction swings slightly, the effects 
on the Doppler beaming of the radiation from this region are dramatic, while balanced by 
the radiation of the larger, slower sheath. This yields a combination of a strongly varying 
component superposed on one having a relatively steady flux density. 

We initially reproduce the i?-band flux density and polarization variations with this 
model, as it is the waveband with the highest sampling during the campaign. We assign to all 
emission elements the average observed optical spectral index of —1.29. To model the degree 
of polarization and flux density of the optical e mission, we utilize t he equations for total flux 
de nsity from opt i cally thin sources deri ved by ICawthornd (120061 ). based on previous work 
by iHughes et al.l (119851 ). The models of ICawthornd (120061 ). created to describe synchrotron 
emission from conically shocked plasma, provide equations for the total flux density from 
regions of both aligned and turbulent magnetic field. The proportionalities, adapted to a 
variable a, are as follows: 



Fiicx^Hsin^) 1 -^- (1) 

and 

Ftu^icnt « Bt~«-(2 - (siiiCb) 1 ~ a (l - * 2 )KCb • (2) 

In the above equations, 5 is the Doppler factor, which is a function of viewing angle 9, 
jet speed /3, and Lorentz factor T: 5 = [r(l — j3sm9)]~ l . The compression factor due to 
shocking in the region is described by factor k and B\\ and -Bturb are, respectively, the parallel 
and turbulent components of the magnetic field. The influence of the magnetic field on flux 
density is described with the angle 9*, the angle between the aligned magnetic field and the 
line of sight, corrected for relativistic aberration. We adapt these formulae for our usage by 
instituting a lack of compression (i.e., no shock, k — 1) and by removing the v factor, as 
we are modeling a single waveband. Applying our observed optical spectral index to these 
formulae, we find the following proportionalities for flux density in the R waveband: 



P ^ r3.29 p>2.29 (n\ 

-TR.spine °spine- D turb,spine > \°J 

-^..transition OC 6\ 3 r ^ tion -Bfj' 29 (sin #*) 2 ' 29 ; (4) 

and 

T? 5:3.29 R 2.29 /r\ 

-PR.sheath ^ "sheath^turb^heath > \° J 

in which we approximate the parallel and turbulent components of the magnetic field to 

be roughly equal for the purpose of simplicity in the model. This approximation is valid 
under the assumption that the magnetic field in the turbulent cells is aligned by shear, 
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but not enhanced. The equations for total flux density from the turbulent sheath and 
spine regions contain factors to account for the strength of the magnetic field and for the 
relativistic beaming effect. The equation for total flux density from the transition region of 
aligned magnetic field contains an additional factor (sin^*) 1- ", reflecting the directionality 
of synchrotron emission, which peaks at a rest frame viewing angle perpendicular to the 
magnetic field. 

Utilizing these equations, we consider a circular cross-section of the jet. See Figure[S]for 
a cartoon of the model structure. For each potential value of the viewing angle (calculated 
for 1000 steps between 0° and 8°), we integrate the flux density across the spine and sheath, 
where the jet cross-section is separated into concentric rings with width 0.1°. We consider a 
range of val ues of the width of the sheath, the width of the spine (within the errors stated 



above from iJorstad et al.ll2005l ). and the average Lorentz factor of the sheath, to find the 
best model fit. We calculate the degree of polarization by approximating that the turbulent 
sheath and chaotic magnetic field component in the spine have negligible polarization, and 
that the boundary region has the maximum polarization possible for synchrotron radiation, 
77.5% for a spectral index a = —1.29. We use this simplifying assumption to minimize 
the number of parameters, given our limited sample size. Further monitoring of this blazar 
would allow the introduction of additional modeling factors, such as the level of magnetic 
field alignment. We then manipulate the viewing angle of the spine and boundary region to 
simulate variations in flux density and polarization. Figure [9] displays the modeled quantities, 
with values discussed below. 

Two factors contribute to the changes in observed flux density with viewing angle. The 
first is the Doppler factor, the decline of which decreases the flux density with increasing 
viewing angle. In addition, the effect on the emission of the ordered component of the 
magnetic field, parallel to the jet axis, varies with the sine factor of equation (4). As the 
de-aberrated viewing angle increases up to 90° - up to an observer's viewing angle in radians 
of ~ 1/r - the observed flux density from the ordered component increases. These two 
multiplicative factors affect the observed emission of the boundary layer such that the flux 
density first grows, then declines, as the viewing angle increases while the spine changes 
direction. The flux density of the turbulent spine is, however, only affected by the Doppler 
factor. The resultant combined flux density for all components peaks at a lower viewing angle 
than does the peak of total polarization. This leaves an observable signature, in which the 
flux density of the combined system can grow to a maximum while the degree of polarization 
rises less steeply and peaks later in time (see Figure E]). 

To model the observed flux density, we find the viewing angle required for an observed 
value of P(-R-band), then determine the modeled flux density corresponding to that viewing 
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angle. We compare the modeled flux density to the observed quantity to find the closest fit 
using our three variables. The best fit is found for a spine of half-width 1.2°, a sheath of 
width 4.6° (between the outer and inner boundaries), and a Lorentz factor in the sheath of 
5.0. Figure M shows the modeled flux density and polarization graphs for this simulation. 
We find that a decrease in viewing angle from 7.5° to 3.2° can create a simultaneous increase 
in flux and polarization, as observed during the 2006 campaign. 

Figure [TD] compares the observed and modeled flux density as well as the observed and 
modeled polarized flux density, and plots the required viewing angle. We demonstrate in 
Figure [lOj panels A and B, the excellent fit of the spine and sheath model to the observed 
flux density and polarized flux density. Following the peak, we observe a decline in flux 
density combined with a delayed decline in percentage polarization. We interpret this as a 
reversal of viewing angle over a period of 7 days from 3.18 ± 0.06° to 7.0 ± 0.4°, as shown in 
Figure [TUl panel C. The optical EVPA (Figure [5]) swings from 6.7° to —37° throughout the 
rise and fall of the flux density. Based on the previously stated relationship between x(core) 
and jet direction, this corresponds to the spine tracing an arc from projected position angle 
—83° to —127° over a 2-week period, with the swing in Xopt providing further evidence for 
rapid, erratic behavior in jet direction. This projected arc swung toward, then away from, 
the line of sight, resulting in the decrease, then increase, of the viewing angle, and thus the 
increase and decrease observed in polarization and flux density. 

We can check the feasibility of such a spine and sheath structure by examining the 
distance downstream necessary for such a boundary layer t o form in the p r esenc e of a velocity 



mismatch. For this purpose, we use a formulation from iJorstad et al.l (120071 ) for shearing 
from transverse velocity gradients. This equation describes the amount of alignment of 
magnetic field expected by such a gradient, assuming that the field is initially turbulent. 
Given the speeds of flow in the spine and sheath, /3 spine and /3 s h ea th, we can calculate the 
fractional alignment f a u g n at a given distance from the initial area of turbulence: 

^ Xf3 le l Az 

/align ~ -p^-, j {&) 

where T is the Lorentz factor of the sheath, /3 rel = (f3 spine - /3 Bhea _ th )(l - ApmeAheath) -1 , R is 
the half-width of the jet, x is the size of the turbulent cells, and Az is the distance between 
the point where the field is completely turbulent and the location of the feature. We set 
an upper limit to R based on component widths close to the 43 GHz core, as determined 
by model fitting of the VLBA data. We estimate x by first calculating the number of cells 
required to reduce the maximum polarization of incoherent synchrotron radiation P ma x to 
the observed value of P within a single reso lution elem ent, N = {P m ax/P) 2 , where N is the 



number of turbulent cells within the beam (jBurnlll966l ). Then x is of order the beam width 
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divided by iV 1 / 3 , ~ 0.049 pc. As a check on our value for x, we calculate the scale size of 
turbulent cells from the 2-day quasi-periodicity in P (optical) shown in Figure El transformed 
to the plasma frame, which gives a value with the same order of magnitude, x ~0.024 pc. 
We find that for the modeled ratio of boundary layer area to turbulent spine area, such 
alignment can occur within 2.5 pc of the core. 

Although the model presented above well describes the observed behavior, we must 
address the issue of a rapidly changing viewing angle, as indicated in Figure [TDJ Previous 
studies of OJ287 have suggested that the viewing angle of the jet is relatively steady over 
short time intervals, with major excursion s only on timescales of years. According to the 
binary black hole model of IValtonen et al.l (120061 ). the viewing angle of the jet has a mean 
value of 2.7° and varies on the order of 0.5° over 30 years. This angul ar variation is too small 
to des cribe the variation of viewing angle in our model. Additionally, iTateyama &: Kingham 
( 120041 ) suggest that the 8 GHz inner jet demonstrates ballistic precession with a period of 
11.6 years. However, this precession is on a time-scale that is orders of magnitude larger 
than the 10 -day variation we m odel. To describe variations in polarized flux and x such as 
we observe, iHughes et al.l (119981 ) consider helical modulation of the velocity vector in OJ287 
by a disturbance in the jet. A helical velocity field would produce the modeled change of 
viewing angle, which swings first toward, then away from, the line of sight. However, due to 
physical constraints, this mechanism requires that rotation of the flow occur on a timescale 
of greater than a year. The variations we report here are too rapid to be explained by either 
a slowly precessing jet or velocity vector rotation. 

Instead, we propose that the change in viewing angle that we infer is most likely the 
combination of wiggles caused by an instability and minor fluctuations in the optical depth 
along the jet. I f the position of the core is defined by t he surface at which the optical depth 
becomes unity (IBlandford fc Konigllll979l ; lKonigllll98ll ). a small change in the optical depth 
of the jet will move the position of the core upstream or downstream. As the core, the 
brightest feature, moves along the jet, it will highlight different parts of the helical velocity 
field. Figure [11] illustrates the jet structure that we propose. If the quasi-steady jet is 
subject to transverse oscillations, perhaps produced by the helical velocity modulation, then 
the change in core position will cause apparent swings in the direction of the section of the 
spine that contains the core. Such movement of the core could be induced by a change 
in magnetic field and electron density in the wake of the disturbance that creates a new 
knot, namely compo nent 2 shown in Figures [lj [2j and [31 This conclusion is supported by 
Jorstad et al.l (120021 ) . in which phase-referencing observations of OJ287 at 43 GHz indicate 
possible motion of the core during a period of prominent flaring. We therefore ascribe the 
swing in viewing angle needed in our model to a temporary change in the position of the core. 
This affects the direction of motion of the emitting plasma in the spine without changing the 
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ly- established long-term integrity of the overall jet. 



previous 

While we are able to use this model to describe the variable emission in 2006 March, 
we must also address the emergence of knots from the core. The polarization vectors of the 
observed knots tend to align well with the direction of their motion. For each component, 
we calculate the difference between its PA and x- We find a difference of 3 ± 9° for the 
combined values of component 2, component 3, and the western component observed only 
in 2005 October. This alig nment can be explained by the propagation of transverse shock 
waves (IHughes et al.l 119851 ). which enhance the magnetic field component in the plane of 
the shock, within the turbulent spine of the inner jet. The presence of relativistic shocks 
within the spine of an otherwise longitudin ally-aligned spine-sheath system has been previ- 
ously theorized (iGopal-Krishna et al.ll2004f) and observed to exist at radio wavelengths (e.g., 
Attridge et al.lll999t iGabuzdal 120031 ) as well as at optical wavelengths on kiloparsec scales 



l e.g.. lDulwich et al.ll2007l ). The observed extent of these jets, as well as the consistently high 
flow speed at large distances, provides evidence for the spine and sheath model. In addition, 
optical observations demonstrate the longevity of magnetic structure in the jet, characterized 
by perpendicular magnetic field in the spine and parallel, aligned magnetic field along the 
boundary layer. 

At times, when the viewing angle and flux density of the shocked emission are such 
that the polarized flux density from the shock exceeds that of the sheared boundary layer, 
the observed x( core ) can switch such that it becomes parallel to the jet. This parallel 
x(core) appeared t hroughout the 1970 's and 1980 's with some rapid switching to and from 
a perpendicular \ Magen-Thornlll980l ). ISillanpaal (119911 ) proposed that changes in \ could 
be due to temporary intervals of dominance of transversely shocke d magnetic fie l d over 
an underlying, structurally stable longitudinal field in the inner jet. iRoberts et al.l (119871 ) 
reported a chang e of y(core) from 160 ° to 99° in late 1982, coincident with the extrapolated 
birth of knot K3 (IGabuzda et al.lll989l) alon g PA = —100°. Our conclusions are supported by 
those of iHagen-Thorn fc Gataullinal (Il99ll ). in which decomposition of the polarized optical 
flux of OJ287 revealed two distinct components, one with \ — 163° (corresponding to our 
sheared component) and the other with x = 84° (corresponding to a transverse shock). 
In our model, the slow and stationary knots observed intermittently in OJ287, such as 
component 1 in Figures [fl [2J and [3] with T = 5.3 ± 0.6c, could be as s ociated with shock s 
propagating along the s lower sheath, as sugge s ted in iKomissarovl (119901 ) , lLaing et al.l (119991 ) , 
Chiaberge et al.l (120001 ). and iGhisellini et al.l (120041 ) . In addition, component 1 exhibits a 
difference between PA and x of 53 ± 15°, suggestive of a more oblique shock than those seen 
in faster components. Thus, the spine and sheath dichotomy can produce distinctly separate 
classes of fast and slow knots as well as the behavior of x observed in OJ287 during our 
campaigns and in previous studies. These include the fast (r = 11.6 to 18.0) and slow (r = 
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0.4 to 1.3) groups of knots noted in OJ287 by IJorstad et aJj (120051 ) 



4.3. Dependence of Degree of Polarization on Frequency and Flux 



We find a clear discrepancy between the optical and radio core polarization in both the 
2005 and 2006 campaigns. In the 2005 observations, P(-R-band) is consistently 2.5 times 
larger - and P(H-band) 2.1 times larger - than P(43 GHz core), while the flux density and 
degrees of polarization are at high levels. In 2006, we observe similar low optical, infrared, 
and radio core percent polarization at a fainter flux density, followed by an increasing ratio 
between the wavelengths as the flux density and P increase (Figure [7]). 

We explain this behavior by appealing to a steepening of the spectrum with frequency. 
We apply a range of spectral indices to the multiwavelength polarization measurements, 
and find a good fit for the following wavelengths: = —1.41, ay = —1.35, = —1.29, 
aj = —1.21, and an = —0.81. These spectral indices reproduce the observed polarization 
ratio changes. The spectral index has little impact on percentage polarization at a high 
viewing angle, which agrees with the similarity in observed values at low flux density. We 
then see a greater impact on polarization at smaller viewing angles, as found in the data, 
with flatter spectra producing a lower percent polarization. To check whether the observed 
steepening in optical index is consistent with synchrotron theory, we examine the indices 
produced by a uniform synchro tron source for a given index of the electron power-law energy 
distribution (jPacholczyki Il970l ). We find that the steepening can be reproduced with a 
synchrotron critical frequency of ~ 1.2 x 10 15 Hz, or A ~ 2500A, for an electron energy 
distribution N(E) oc E~ s with index s = 2.5. 

We examine the optical spectra to determine whether this steepening in index was in 
fact observed. We are not able to detect such steepening within the optical region given 
the error in spectral index for individual wavebands, cr(o;)~ ±0.25. In Figure [121, we 
display the average optical-to-infrared spectrum overlaid with the model steepened spec- 
trum. For each epoch at which both infrared and optical measurements are available, 
we normalize the observed spectrum to log [Flux Density(ff)]=l, then average the spec- 
tra to produce the displayed data. As is shown in Figure [12l we observe a clear steep- 
ening between the infrared to far-red index = —0.9 ± 0.1) and the optical index 
[ot opt = —1.3 ± 0.1), which is fit well by the model steepened spectrum. In addition, all 
five modeled indices agree with the observed values within one standard deviation. Histori- 
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When we extend the spectrum to fit the 43 GHz data, we find that the percentage 
polarization in the radio requires an optically thin spectral index near zero. From this, 
we conclude that the 4 3 GHz po l arizat ion is diluted by emission from an optically thick 
region. As described in iGabuzdal (120031 ). the agreement between x(optical) and x(43 GHz 
core) suggests that the polarized radio emission is dominated by the optically thin region. We 
calculate t he necessary dilutio n by the weakly polarized optically thick zone using the formula 
derived in iPacholczyk! (jl970l ). under the assumption that the optically-thin spectral index 
~ —0.8, ro ughly that inferred fo r the i nfrared iJ-band (cf. the spectral energy distribution 
derived by llmpey fc Neugebauerl (119881 )). From the observed percentage, we determine that 
~ 30% of the 43 GHz core emission is optically thick. 



4.4. Polarization Behavior During 2005 Campaign 

Examination of the 2005 campaign results reveals behavior that, while less dynamic, 
agrees with the above model. As Figure H] indicates, we observe an upward trend in both 
optical and 43 GHz core flux, as well as an increase in P(optical), accompanied by a swing 
in x(optical). A change in viewing angle of the spine can explain all three observed effects. 



4.5. Alternative Explanations of the Polarization and Variability 



We consider two alternative expl anations for the observ ed behavior of \ in the 43 GHz 
core of O J287. In one recent proposal, iGabuzda et al.1 (120041 ) use complex rotation measures 
and close examination of % values to describe the parsec-scale jets of BL Lac objects as 
containing a helical field in which the toroidal component of magnetic field dominates. While 
this description is consistent with the steadiness of x i n components 2 and 3 as they proceed 
down the jet, it fails to explain either the dynamical nature of the core or the fast increase 
in flux and polarization. I n addition, although many BL Lac objec ts display a %(core) 
parallel to the jet direction (IGabuzda et al.ll2000l ; iMarscher et al.ll2002l ). our study of OJ287 
shows a highly-polarized x(core) perpendicular to the jet, a condition that would necessitate 
extreme, nearly longitudinal pitch angles of the helix. Historically, the core of OJ287 has 
high polarization in instances when it maintains a transverse y - shown in our data - as we ll 
as when it possesses a longitudinal x (|Hagen-Thornlll980l ; lHagen- Thorn &: Gataullinalll99ll ). 
which is not expected in the helix model. 



Another study of OJ287, by IVicente et al.1 (119961 ) using the models of iHarded (119871 ). 
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suggests that the progression of components along the jet is in fact the motion of shocks 
following a helical path. If this path were to change by 90° within the resolution of the 
VLBA, it could very well describe the dichotomy we see between \ in the core and that of 
its ejected components, as well as the change in viewing angle as plasma proceeds down the 
parsec-scale jet. Nonetheless, we observe no helical behavior in the path of motion of the 
components within pos i tiona l errors (see Figure [3]), contradicting the model. Furthermore, 



Tateyama fc Kinghaml (120041 ) have examined this model for OJ287 and conclude that it is 



insufficient to fully describe the data, proposing instead that ballistic precession produces a 
changing inner jet direction on lo nger timescale s than seen in our observations. We note, 
however, that the conclusions of IVicente et all (119961 ) described components on scales of 
several milliarcseconds and on timescales of several years, and are therefore not entirely 
applicable to the data that we display. 



5. Summary and Conclusions 

The increase of polarization across millimeter to optical wavelengths that our observa- 
tions reveal in OJ287 leads us to conclude that the optical emission is cospatial with the 
43 GHz core. This conclusion is supported by the similarity in value and behavior with 
time of the EVPA at different wavebands. The EVPA of the core has maintained an angle 
roughly perpendicular to the flow of the jet both historically - except for occasional 90° flips 
- and in our observations. This steadiness implies persistence of a region of velocity shear 
in the core region. We have proposed a model in which velocity shear creates an ordered 
longitudinal component of the magnetic field in the transition region between a fast, turbu- 
lent spine of the jet, and a slower, turbulent sheath. Shock waves emerge along the spine 
to form superluminal knots, and along the sheath to produce relatively slow, or stationary, 
features. We model the core with this spine-sheath structure, with the inner jet containing 
flux concentrations. With the existence of velocity vector gradients in the inner jet, changes 
in the position of the core cause variations in position and viewing angles. The occasional 
appearance of shocks, whose emission dominates the flux density, can explain the 90° jumps 
in x(43 GHz core) observed at previous epochs. 

Our observations of OJ287 have revealed a dynamic jet downstream of the core. It 
no longer shows relatively slow and steady motions as observed in the 1980's and 1990's, 
nor does it display non-core stationary features as observed in the late 1990's. We have 
uncovered a tendency for the jet to produce knots following different trajectories at intervals 
of months. The viewing angle changes that we propose are counter to those called for in 
studies of slower, systematic viewing angle changes in OJ287. Rather, they correspond to 
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more rapid fluctuations in the direction of the flow of the emitting plasma within the core. 
We speculate that disturbances propagating down the jet cause shifts in the location of the 
site where the optical depth equals unity, which is identified as the core on VLBA images. If 
the direction of the spine is subject to an instability that results in spatial oscillations, the 
shifting core position would appear as a swing in the direction of the fast spine. This change 
in angle between the flow velocity vector and the line of sight would then be responsible for 
the variability in flux and polarization that we have observed. 

We have collected similar data from a number of other blazars during the 2005 October- 
November campaign and the 2006 March-April campaign. We are currently analyzing the 
data to determine the level of generality of the conclusion that polarized optical emission is 
correlated in variability and polarization properties with the 43 GHz core in blazars. The 
result of this analysis, in concert with previously established correlations between optical 
emission and high-frequency emission, will contribute toward the development of a general 
model for emission in parsec-scale jets. This will be described in a summary paper to follow. 
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Epoch 


Component 


R (mas) 


PA (°) 


1 (Jy) 


P (%) 


X° 


L)\J r 


05 Oct. 24 


Core 






2.01 


12.7 


21 


13.7 




West Component 


0.08 


-58 


1.09 


7.0 


-68 


5.2 




Component 1 


0.32 


-147 


0.03 


0.0 




3.1 


05 Oct. 28 


Core 






2.30 


11.8 


16 


9.6 




West Component 


0.09 


-84 


0.77 


6.1 


-87 


3.1 




Component 1 


0.38 


-154 


0.02 


0.0 




2.0 


05 Nov. 2 


Core 






2.33 


11.8 


19 


9.1 




West Component 


0.11 


-80 


0.65 


7.6 


-83 


3.9 




Component 1 


0.40 


-159 


0.03 


0.0 




2.6 


06 Mar. 27 


Core 






1.55 


9.7 


-4 


41.8 




Component 2 


0.10 


-90 


1.35 


0.0 




28.2 




Component 1 


0.54 


-145 


0.03 


69.9 


-8 


20.4 


06 Mar. 31 


Core 
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13.8 


1 


15.7 




Component 2 


0.12 


-91 


2.01 


0.8 


75 


4.6 




Component 1 


0.60 


-145 


0.03 


61.2 


-22 


4.2 


06 Apr. 5 


Core 






1.82 


14.4 
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1 Chi 2 ~ )OF = Chi 2 j (2 x no. visibilities - no. modclfitting parameters) of the model containing the respective 
component plus the components above it. 
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Fig. 1. — 43 GHz VLBA images of OJ287 at 6 epochs of observation (dates listed). The bars 
represent polarized intensity, and the position angle of the bar corresponds to the EVPA. 
Neither Faraday correction of the core EVPA nor statistical bias correction of the polarized 
intensity has been applied to this display of the images. The contour levels are 0.5, 1, 2, 
4, 8, 16, 32, and 64% of the peak total intensity. Left: The restoring beam (shown in the 
lower right corner) has dimensions of 0.33 mas x 0.16 mas at an angle of —18.7°. Right: The 
high-resolution restoring beam (shown in the lower right corner) has a circular diameter of 
0.1 mas, and is representative of the resolution of the longest baseline. 
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Fig. 2. — 43 GHz VLBA images of OJ287 at 3 epochs of observation for (Top) the 2005 
campaign and (Bottom) the 2006 campaign (dates listed). The contour levels are 0.5, 1, 2, 
4, 8, 16, 32, and 64% of the peak total intensity. The restoring beam (shown in the lower 
right corner) has a circular diameter of 0.1 mas, and is representative of the resolution of 
the longest baseline. 
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Fig. 3. — Left: Separation in milliarcseconds relative to the core of three components, with 
one point plotted for each epoch at which the component appears on the respective VLBA 
image. The triangles (1) correspond to the component with birth date ~2004 March, the 
filled circles (2) to the component with birth date ~2006 February, and the open squares (3) 
to the component with birth date ~2006 May. Robust least-squares fitting to the separation 
of each component is shown with solid lines. Julian date 2453700.5 corresponds to 2005 
November 26. Right: Map of the locations of the three components. Robust least-squares 
fitting to the position angle of each component is shown with solid lines. The motion is 
outward from the core, located at position (0,0). 
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Fig. 4. — Data from the 2005 campaign for O J287. Left: Data in R-b&n d as a function of time 
showing flux density (A) corrected for reddening (ISchlegel et al.lll998l ). percent polarization 
(B), and \ (C)- The dashed line represents the average x(43 GHz core). Right: Data for 
the 43 GHz core displayed in a similar manner as the optical data. In all cases, x( core ) has 
been corrected by 30° to account for Faraday rotation. The solid line represents the average 
x(R)- Julian date 2453668.5 corresponds to 2005 October 25. 
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Fig. 5. — Measurements from the 2006 campaign for OJ287 shown in the same manner as in 
Figure 4, except for the addition of if -band measurements (in (A), (B), and (C)) represented 
by a dash-dotted line. Julian date 2453820.5 corresponds to 2006 March 26. 
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Fig. 6. — Multiwavelength plot of polarization variations with time during the 2005 cam- 
paign. B waveband represented by dotted line, V by dashed line, R by solid line, / by 
dash-triple dotted line, and H by dash-dotted line. Julian date 2453668.5 corresponds to 
2005 October 25. 
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Fig. 7. — Multiwavelength plot of polarization variations during the 2006 campaign. Julian 
date 2453820.5 corresponds to 2006 March 26. 
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SPINE: 




•Turbulent 






•Can contain shocks 







•P = 0% 






•r = 16.5 






•Half-width = 1 .2° 







SHEATH: 




•Turbulent 




•Can contain shocks 


• 


•P = 0% 




•r = 5.0 




•8(outer)-9(inner) = 4.6° 



BOUNDARY: 
•Shear aligned 
•P = 77.5% 
•r = 10.75 

•0(outer)-9(inner) = 1.4 



Fig. 8. — Cross-sectional diagram of modeled jet in the vicinity of the core, with nominal 
quantities for the spine, sheath, and boundary regions listed. 
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Fig. 9. — Modeled flux density (top) and degree of polarization (bottom) as a function of 
viewing angle. Flux density has been normalized to the peak. Polarization peaks at a viewing 
angle of 4.05°, and flux density at 0.0°. 
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Fig. 10. — Model fit to 2006 local flux density peak. A: Observed flux density at R waveband 
(filled circles) and modeled flux density (open diamonds). B: Observed polarized flux at R 
waveband (filled circles) and modeled polarized flux density (open diamonds). C: Modeled 
viewing angle of the jet in degrees. Julian date 2453820.5 corresponds to 2006 March 26. 
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Fig. 11. — Oblique view of jet. Core position shifts back and forth, causing swings in the 
velocity vector of plasma in the core. Line segments represent direction of the magnetic field, 
which is parallel to the velocity vector in the boundary layer. Viewing angle of spine and 
boundary layer vary with time. 
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Fig. 12. — Optical-to-infrared spectrum from 2006 campaign. Displayed above is the average 
spectrum, normalized to peak at unity in the if -band. The optical spectrum is binned in 4 
equal segments of wavelength between 4200A and 7500A. The solid line displays the modeled 
spectrum. 



